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A simple approach for computing acceleration and velocity of a structure from the strain 
is proposed in this study. First, deflection and slope of the structure are computed from the 
strain using a two-step theory. Frequencies of the structure are computed from the time 
histories of strain using a parameter estimation technique together with an autoregressive 
moving average model. From deflection, slope, and frequencies of the structure, acceleration 
and velocity of the structure can be obtained using the proposed approach. Simple harmonic 
motion is assumed for the acceleration computations, and the central difference equation with 
a linear autoregressive model is used for the computations of velocity. A cantilevered 
rectangular wing model is used to validate the simple approach. Quality of the computed 
deflection, acceleration, and velocity values are independent of the number of fibers. The 
central difference equation with a linear autoregressive model proposed in this study follows 
the target response with reasonable accuracy. Therefore, the handicap of the backward 
difference equation, phase shift, is successfully overcome. 


Nomenclature 

AR = autoregressive 

ARMA = autoregressive moving average 

DOF = degrees of freedom 

FBG = fiber Bragg grating 

FE = finite element 

FOSS = fiber optic strain sensor 

NASA = National Aeronautics and Space Administration 

SEREP = System Equivalent Reduction and Expansion Process 

a i i, a 2 i = autoregressive coefficients for the i-th orthonormalized coordinate rj t 

k = discrete time 

At = time step 

r]i(k) = i-th orthonormalized coordinate at discrete time k 

cOi = i-th frequency from measured strain time histories 

{ Q a )k = generalized aerodynamic load vector at discrete time k 

{q} k = generalized coordinates vector at discrete time k 

{ qjyi}k = master degrees of freedom vector at discrete time k 

{q M } k = measured master degrees of freedom vector at discrete time k 

{ q s } k = slave degrees of freedom vector at discrete time k 

{e} k = strain vector at discrete time k 

{q} k = orthonormalized coordinates vector at discrete time k 

[G] = damping matrix 

[K] = stiffness matrix 

[M] = mass matrix 

[O] = eigen-matrix 

[*m] = eigen-matrix corresponds to master degrees of freedom 

[<*>*] = eigen-matrix corresponds to slave degrees of freedom 

{*} = the first time derivative of vector {*} 
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{*} = the second time derivative of vector {*} 

[*] _1 = inverse of a matrix [*] 

[*] r = transpose of a matrix [*] 

I. Introduction 

F LEXIBLE and light weight optical fibers not only revolutionize the telecommunications, but also innovate the 
sensing world. Optical fibers can be used as fiber optic sensors to measure strain, temperature, pressure, and other 
quantities. Properties of light, such as intensity, phase, transit time of light, et cetera, that propagate in the fiber can be 
altered by modifying a fiber. Recently, fiber optic sensors have been developed to measure co-located strain 
simultaneously with very high accuracy using fiber Bragg gratings (FBGs). Specifically, a fiber optic strain sensor 
(FOSS) uses a series of FBGs to obtain measurements at intervals as small as every millimeter 1 along a fiber and a 
frequencies of several kHz. 2 The ability of FBGs to operate at such high frequencies makes them an ideal choice for 
both static and dynamic aerospace applications. 

Obtaining wing deflections during flight are essential measurements for structural health monitoring and active 
flexible motion controls. Calculating wing deflection values from the measured FOSS data is available along the 
singular optical fibers. Ko et al. 3 uses closed-form equations to calculate real-time operational loads and deflections 
on complex structures using surface strain data. The Kang et al. 4 calculation method uses a matrix equation based on 
the displacement and strain mode shape matrices. The method described by Moore 5 includes the steps of deducing 
curvature and bend radius data for the optical fiber from the strain data: curve-fitting to obtain curvature and bending 
direction functions pertaining to the optical fiber, calculating a torsion function using the bending direction function, 
and calculating the three-dimensional shape of the optical fiber from the combination of the curvature, and bending 
and torsion functions. These three methods are suitable for calculating deflection and load values in cross-sectional 
lumps based on the location of the sensors used to measure strain on the surface. 

The availability of wing deformation at all element grid points across the structural finite element (FE) model will 
improve the ability to monitor the health and behavior of a complex three dimensional structure such as an aircraft in 
flight. Several studies have been performed which evaluated using an inverse finite element method. 6,7 ’ 8 This method 
is based on the minimization of a weighted-least-squares functional that uses the complete set of strain measurements. 
The two-step approach developed by Pak 9 computes wing deflection values along the location of all the FOSS first, 
and then these deflection values are expanded to the entire structure using the System Equivalent Reduction and 
Expansion Process (SEREP) 10 technique. 

Acceleration and velocity of structural motion are also needed to compute internal and external loads during flight 
through the use of the structural governing equations of motion. This paper focuses on a computation of acceleration 
and velocity over the entire three dimensional structure based on measured strain information. First, unsteady 
structural deformation of the entire three dimensional structure is obtained using the two-step approach introduced by 
Pak 9 . Next, accelerations and velocities are computed using a simple harmonic motion assumption and a central 
difference equation, respectively. 


II. Mathematical Background 

A long term objective of the current study is the computation of external aerodynamic loads and internal structural 
loads from strains which are measured from a FOSS. This objective can be accomplished through the use of the 
structural dynamic governing equations of motion. Mass, damping, and stiffness matrices in Eq. (1) are from a finite 
element model. From Eq. (1), deflection, acceleration, and velocity values are needed to compute elastic, inertia, and 
damping loads. Therefore, computing deflection, acceleration, and velocity information from the measured strain are 
the critical quantities for the long term objectives. In the first section, deflections and slopes of an entire structure are 
calculated from measured strain through the use of a recently developed two-step theory. 9 Acceleration and velocities 
are computed in the second and third sections using a harmonic motion assumption and the central difference equation, 
respectively. 

A. Computation of Wing Deflection from Measured Strain 

Consider the following structural dynamic governing equations of motion shown in Eq. (1): 

[M m k + [G]{q} k + [K ]{q} k = (1) 

where [M], [G], and [K] are mass, damping, and stiffness matrices, respectively, and {q} k and {Q a }k are the 
generalized coordinates and the generalized aerodynamic load vectors at discrete time k , respectively. 
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Deflections along the location of the FOSS can be computed from the measured strain data { e} k using a piecewise 
least-squares method, an Akima spline, and a linear assumption as described in the two-step approach. 9 These 
computed deflections at discrete the FOSS locations are combined with a FE model of the structure in order to 
interpolate and extrapolate the deflection and slope of the entire structure through the use of SEREP. Re-arranging all 
the degrees of freedom (DOF) in the FE model yields Eq. (2): 


iqh = [ q q M s } k = [*]fa} k 




ivh 


( 2 ) 


where {q M } fc is the master DOF at discrete time k. In this approach, deflections along the fiber locations computed 
from the first step of the two-step approach 9 are defined as the master DOF. The rest of the deflections and slopes all 
over the structure are defined as slave DOF at discrete time k , {qs)k- In Eq. (2), matrices [d> M ] and [0 5 ] are 
eigen-matrices corresponding to master and slave degrees of freedom, respectively, and [q] k is the orthonormalized 
coordinates vector at discrete time k. Therefore, the following two equations, Eqs. (3) and (4), are derived from Eq. (2). 


— MW* P) 

{Qs}k = [<t>s]{ri}k (4) 

In Eq. (3), changing master DOF at discrete time k , to the corresponding measured values gives Eq. (5): 

{ qjyi)k — ( 5 ) 

where {q M } fc is obtained from the first step of the two-step approach. 9 Pre-multiply [0 M ] r to Eq. (5) with matrix 
inversion gives Eq. (6) for computing the orthonormalized coordinates vector at discrete time k, 

[qh = {[^MVi^r^MViqMh, (6) 


and the generalized coordinates vector {q} k is obtained from Eq. (2). 


B. Computation of Acceleration from Computed Wing Deflection and Frequencies 

From Eq. (2), acceleration at discrete time k can be expressed as shown in Eq. (7). 



O 
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{r)}k 


(7) 


By assuming simple harmonic motion for the i-th orthonormalized coordinate, rji(k) = —ojfrjiik) i = 1,2, 
acceleration can be obtained from deflections and frequencies as shown in Eq. (8). 
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Substituting Eq. (8) into Eq. (7) gives Eq. (9). 




’m [<**?] 

sNI 


( 8 ) 


(9) 


Therefore, acceleration is computed using Eqs. (6) and (9). In Eq. (8), frequencies co* i = 1,2 , ... ,n , are computed 
from unsteady strain distribution at a selected point using an on-line parameter estimation technique together with an 
autoregressive moving average (ARMA) model. 11 
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C. Computation of Velocity from Computed Wing Deflection and Frequencies 

From Eq. (2), velocity at time k can be written as shown in Eq. (10). 



( 10 ) 


Consider the following central difference equation, Eq. (11), for the time derivative of orthonormalized coordinates 
vector {i]} k . 


{j>}fc = 


fo } fc + 1 - Mfc - 1 
2 At 


The time derivative of the i-th orthonormalized coordinate can be written in Eq. (12). 


(ii) 


MW 


ViXH + 1) ~ ViXk - 1) 

2At 


( 12 ) 


It should be noted in Eqs. (11) and (12) that the orthonormalized coordinates vector at future time k+\ is needed to 
have the time derivative of the orthonormalized coordinates vector at time k. 

Consider a linear autoregressive (AR) model 11 for the i-th orthonormalized coordinate 77 * (fc) given in Eq.(13): 

Vi(.W = a^riitk - 1) + a 2i rh(k - 2) (13) 


where AR coefficients for the i-th orthonormalized coordinate, a lt and a 2 u are computed from the same i-th frequency 
co i in Eq.(8). From Eq. (13), the i-th orthonormalized coordinate at future time k+\ can be expressed as shown in 
Eq.(14). 


V i(k + 1) = a^Thik) + a 2 jJ7i(/c - 1) 


Substituting Eq. (14) into Eq. (12) gives Eq.(15). 


MW 


OliViOO + ( a 2i ~ l)?7i(fc ~ 1) 

2At 


In summary, velocity of a flexible wing motion is computed as follows: 


(14) 


(15) 


Step 1: compute the orthonormalized coordinates vector {t]} k from Eq. (6). 

Step 2: compute the time derivative of the orthonormalized coordinates vector {i]} k from Eqs. (15) and (16). 

piiW) 

{«» = ffH (16) 

Urn) 


Step 3: compute the time derivative of the generalized coordinates vector { q} k from Eq. (10). 


Therefore, deflection, acceleration, and velocity based on the generalized coordinates vectors { q } k , { q } k , and { q } k , as 
well as the orthonormalized coordinates vectors {t]} k , rj } k , and {rj} k , are computed from the measured strain data 

to*- 
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III. Results and Discussions 

The theory presented in this study is validated using a cantilevered rectangular wing model as shown in Fig. 1. 
This FE model represents a wind tunnel test article 12 that has a 0.065-inch thick aluminum plate with a thin plastic 
foam skin that forms a 6% circular arc cross-sectional shape. A transient load using a half sine short duration pulse is 
applied at the tip leading edge corner. Strain values are output from each zero stiffness beam element segment at 
simulated FOSS fiber locations along the span, a total of 50 beam elements for each fiber and along upper and lower 
plate elements as shown in Fig. 2. The upper and lower elements are offset from the center plate elements to generate 
non-zero strain. Note that these elements have a finer distribution at the root as required to lower the propagation of 
integration errors outwards to the tip. 

For validation and refinement of this theory, a computational validation was performed. This wing model is 
analyzed using the MSC/NASTRAN. 13 The computed strain values are assumed to be strain output from a FOSS. The 
computed deflection, acceleration, and velocity values are considered as the targets to be compared with the results 
obtained from the current approach. 

The cantilevered rectangular wing FE model has two versions that differ in material and mass properties to test the 
current approach with a real world variation of model accuracy as compared to measured frequency results. These two 
FE model versions are shown in Table 1 where case 1 has its plastic foam mass equally smeared among the aluminum 
plate density, and case 2 has the plastic foam mass distributed as lumped masses over the 6% circular-arc cross 
sectional shape. Young’s modulus and shear modulus for the case 2 model are tuned using the structural dynamic 
model tuning tool developed at the National Aeronautics and Space Administration (NASA) Armstrong Flight 
Research Center. 14 In Table 2, computed modal frequencies for cases 1 and 2 are compared to available measured 
values. 12 The case 1 frequencies have up to a 5.6 % error which is typical of such FE models near the limits of the 5% 
frequency NASA FE modeling standard. 15 The frequency errors for case 2, at no more than 0.8%, represent an FE 
model that has the best accuracy that can usually be expected. 

Both cases 1 and 2 were run to generate the target responses with NASTRAN modal transient response analysis 
(sol 112) with 1200 time steps. The case 1 mode shapes are used to calculate transformation matrices for both cases. 
The case 1 mode shapes are the eigen and comparison functions for the case 1 and case 2 transient responses, 
respectively. 

Frequency estimations are generated using computed strain values for cases 1 and 2 and are taken at the leading 
edge root as shown in Fig. 2. The computed strain time histories for both cases are shown in Fig. 3. The frequency 
estimation uses Bierman’s U-D factorization algorithm with 20 AR coefficients, a covariance matrix resetting interval 
of 80 steps, and a forgetting factor of 0.98. For cases 1 and 2, sampling times of 0.0006267 s and 0.0006487 s, and 
Nyquist frequencies of 797.9 Hz and 770.8 Hz, respectively, were used. These estimated frequencies are compared to 
the target values generated from the NASTRAN output shown in Table 3. The estimated frequency errors are very 
low, if not at zero, with the highest standing out at -0.09% for case 2, mode 1. 

Using the current approach, deflections and pitch angles are calculated using computed strain along simulated fiber 
locations of the wing model. For case 1, target deflection, acceleration, and velocity values are compared to the current 
approach calculations using 22 fiber locations at the wing model leading edge tip grid 51, shown in Fig. 2. Since grid 
51 is located on the leading edge fiber, its deflections are directly computed from strain values and are master DOF. 
Using wider and close up time histories, the out of plane, or z deflections are shown in Fig. 4, and the pitch angles are 
shown in Fig. 5. The green dot line and the black solid line represent the target and values from the current method, 
respectively. The z deflection comparisons are so close to each other that no discernable differences can be seen in 
the plots between the target and current method values. The pitch angle comparisons are almost as equally 
indiscernible from each other except at several of the areas with large changes in gradient where, in general, the target 
values are a little larger than those for the current method. 

In Fig. 4(a), the period of the large oscillation is approximately 0.065 s, 15.4 Hz. Therefore, the first natural 
frequency is observed in the z deflection response. On the other hand, seven oscillations are observed for 0.09 s in 
Fig. 5(a). Therefore, the period of the large oscillation is approximately 0.0129 s, 77.8 Hz. This frequency corresponds 
to the second mode, 77.4 Hz, and is the first torsion mode. In Fig. 5(b), the small oscillation is also dominant, and 
three oscillations are observed during 0.0041 s, and therefore the corresponding frequency is approximately 730 Hz, 
period is 0.00137 s. 

Out of plane, or z direction, accelerations are compared in Fig. 6, and comparisons for pitch accelerations are 
shown in Fig. 7. Both acceleration plots are showing an initial opposite trend in phase between the target and current 
method values for about a half cycle, but the current method very quickly recovers and begins to match up closely. 
This initial phase difference is mainly due to the simple harmonic assumption for the transient responses. For the z 
accelerations, the comparisons are showing several peak values with slightly larger target values over those for the 
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current method. The pitch acceleration comparisons are showing the same match of trend with more of the target peak 
values showing some larger values over those for the current method. 

Out of plane velocity and pitch rate based on the current method values are compared with the corresponding target 
responses in Figs. 8 and 9, respectively. Both figures show that the current method values are smaller than the target 
values during the transient time period. Slightly larger target values over the corresponding values for the current 
method are also observed in the z velocity and pitch rate responses. 

The case 2 wing model comparisons of the target results to the current method using 22 fiber locations are made 
at a tip 1/3 aft chord location, grid point 2601 as shown in Fig. 2. The z deflections and pitch angle time histories are 
shown in Figs. 10 and 11, respectively. The matching and small deviations of the target values to those of the current 
method are essentially the same as they were for case 1, for both the z deflections and the pitch angles. In Fig. 10(a), 
the period of the large oscillation is approximately 0.068 s, 14.7 Hz. Similar to case 1, the first natural frequency 
14.29 Hz is observed in the z deflection response. On the other hand, eight oscillations are observed for 0.10 s in 
Fig. 1 1(a), at approximately 80 Hz. Therefore, the 80.18 Hz torsion mode is dominant in pitch motion. 

The z and pitch accelerations time history comparisons are shown in Figs. 12 and 13, respectively. The z velocity 
and the pitch rate comparisons are given in Figs. 14 and 15, respectively. And again, these results are essentially 
repeating the same trends as seen with case 1 . 

The target to the current approach results for the wing model cases 1 and 2 are summarized in Table 4. An additional 
sensitivity study of a fewer number of fibers for case 2 is also shown in Table 4. The previous number of fibers used 
for case 2 is lowered from 22 to 10 and 6, and is depicted in Fig. 16. Percent errors between the target and current 
method results are calculated using Eq. (17). 


% Error = 


Yk=o \ Current approach (fc) - Target(/c)| 

2fc=o I T ar g et (fc) I 


(17) 


The target values are those directly computed from the FE wing model. The results at discrete time, k are summed for 
all the time histories. 

Table 4 is showing a very small general increase in error with a reduction of the number of fibers for case 2. This 
minimal increase in error is an indication of the robustness of the least-squares surface fitting procedure from the fiber 
locations out to the rest of the wing model. For case 1, with the z and pitch accelerations at grid 51 and on a fiber 
location, the errors in z acceleration are noticeably lower than those at an off fiber location at grid 2601. It should be 
noted that the z deflection at grid point 51 belongs to the master DOF. However, any DOF at grid point 2601 belong 
to the slave DOF. The z deflection results have the lowest percent error, at no more than 1.79%. The pitch angles, z 
accelerations, pitch accelerations, z velocity, and pitch rate calculated during and after the second step, as slave DOF; 
have higher errors are higher at 6.35%, 17.6%, 10.2%, 19.0%, and 12.0 %, respectively. It should also be emphasized 
that the comparison functions are used for case 2 instead of using eigen functions. However, the quality of the results 
are not degraded too much. 

Finally, instead of using the central difference equation, Eq. (11), the following backward difference equation 
given in Eq. (18) is used for the pitch rate computations with the case 2 model. 


{j)}k = 


0?}k - Wk-i 

At 


(18) 


Comparisons between the central difference equation with a linear AR model and the backward difference equation 
are shown in Fig. 17. In Fig. 17(b), small phase differences are observed between the target responses and the values 
from the current study. It should be concluded that the current proposed central difference equation with a linear AR 
model closely follows the target response. 


IV. Conclusion 

A simple approach for computing acceleration and velocity of a structure from the time histories of strain is 
proposed in this study. First, deflection and slope of the structure are computed from strain using a two-step theory. 
Next, the simple harmonic motion assumption is used for the acceleration computations. Therefore, accelerations are 
computed from the deflections obtained from the first step of the two-step theory and the estimated frequencies. 
Frequencies of the structure are computed from the time histories of strain using a parameter estimation technique 
together with an autoregressive moving average model. Finally, the central difference equation with a linear AR model 
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is used for the computations of velocity. A cantilevered rectangular wing model is used to validate the simple 
approach. 

Quality of the computed deflection, acceleration, and velocity values are independent of the number of fibers 
(number of master DOF) except using too small of a number. If only a single fiber is used for the whole surface, then 
the torsional motion is not captured properly. 

Velocity values based on the central difference equation with a linear AR model proposed in this study follows the 
target response with reasonable accuracy. Therefore, the phase shift issue associated with the backward difference 
equation is successfully overcome using the central difference equation with a linear AR model proposed in this study. 

A future application of the approach will be performed on the X-56A aeroelastic research aircraft (Lockheed 
Martin, Bethesda, Maryland) in order to generate pre-flight test estimates of aerodynamic forces on the wing surfaces. 
These force estimates will be correlated to those results as derived from onboard FOSS strain flight data. 
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Tables 

Table 1. Wing finite element model properties. 


Properties 

Case 1 model 

Case 2 model 

Young’s modulus 

9847900 psi 

9207766 psi (tuned value) 

Shear modulus 

3639672 psi 

3836570 psi(tuned value) 

Density 

0.11166 lb/inch 3 

0.1 lb/inch 3 

Foam weight 

Smeared 

Lumped mass 

Total weight 

0.3806 lb 

0.3806 lb 

Xcg 

2.28 inch 

2.28 inch 

Yog 

5.75 inch 

5.75 inch 

Thickness 

0.065 inch 

0.065 inch 


Table 2. Measured versus computed frequencies. 


Mode 

Measured 12 (Hz) 

Case 1 (Hz) 

% Error 

Case 2 (Hz) 

% Error 

1 

14.29 

15.09 

5.6 

14.29 

0.0 

2 

80.41 

77.40 

-3.7 

80.17 

-0.3 

3 

89.80 

93.57 

4.2 

89.04 

-0.8 

4 

N/A 

246.37 

N/A 

248.76 

N/A 

5 

N/A 

262.02 

N/A 

252.41 

N/A 

6 

N/A 

455.22 

N/A 

459.34 

N/A 

7 

N/A 

511.27 

N/A 

485.61 

N/A 

8 

N/A 

642.72 

N/A 

606.65 

N/A 

9 

N/A 

722.32 

N/A 

718.59 

N/A 

10 

N/A 

773.93 

N/A 

747.65 

N/A 


Table 3. Estimated system frequencies for cases 1 and 2. 


Mode 

Case 1 

Case 2 

Target (Hz) 

Estimated (Hz) 

% Error 

Target (Hz) 

Estimated (Hz) 

% Error 

1 

15.09 

15.09 

0.00 

14.29 

14.28 

-0.09 

2 

77.40 

77.40 

0.00 

80.17 

80.18 

0.02 

3 

93.57 

93.57 

0.00 

89.04 

89.05 

0.01 

4 

246.37 

246.37 

0.00 

248.76 

248.77 

0.00 

5 

262.02 

262.02 

0.00 

252.41 

252.41 

0.00 

6 

455.22 

455.22 

0.00 

459.34 

459.34 

0.00 

7 

511.27 

511.27 

0.00 

485.61 

485.61 

0.00 

8 

642.72 

642.72 

0.00 

606.65 

606.65 

0.00 

9 

722.32 

722.32 

0.00 

718.59 

718.60 

0.00 

10 

773.93 

773.93 

0.00 

747.65 

747.66 

0.00 


Table 4. Target and current method errors for case 1 and 2 simulations. 


Model 

Grid (# of fiber) 

% Error 

Deflection 

Acceleration 

Velocity 

Z 

Pitch 

Z 

Pitch 

Z 

Pitch 

Case 1 

51 (22) 

1.55 

5.36 

6.42 

7.96 

10.5 

12.0 

Case 2 

2601 (22) 

1.38 

5.76 

16.9 

9.84 

15.0 

11.4 

2601 (10) 

1.67 

5.99 

17.0 

10.2 

15.9 

11.7 

2601 (6) 

1.79 

6.35 

17.6 

10.2 

19.0 

11.8 
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Figures 

22 Simulated FOSS locations using beam elements 
Strain data is taken at each of the 50 beam elements per fiber. 



Figure 1. Cantilevered rectangular wing model. 



Figure 2. Wing model with simulated FOSS locations. 
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a) Case 1 . 


b) Case 2. 


Figure 3. Strain time history at the leading edge of the wing root section for case 1 and case 2. 


0.3 



0.00 0.02 0.04 0.06 0.08 0.10 

Time (sec) 


a) From 0 to 0.10 s. 


0.3 

: Target 

0.2 : Current Method 


t o.i 



0.00 0.01 0.02 0.03 

Time (sec) 

b) Zoom in from 0 to 0.03 s. 


Figure 4. Case 1, Z deflection time history at grid 51. 
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Figure 5. Case 1, pitch angle time history at grid 51. 
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Figure 6. Case 1, Z acceleration time history at grid 51. 
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Figure 7. Case 1, pitch acceleration time history at grid 51. 
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Figure 8. Case 1, Z velocity time history at grid 51. 
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Figure 9. Case 1, pitch rate time history at grid 51 
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Figure 10. Case 2, Z deflection time history at grid 2601. 



- 0.05 


0.00 0.02 0.04 0.06 

Time (sec) 

a) From 0 to 0.10 s. 


0.08 


0.05 

0.04 


0.10 



Time (sec) 

b) Zoom in from 0 to 0.03 s. 


Figure 11. Case 2, pitch angle time history at grid 2601. 
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Figure 12. Case 2, Z acceleration time history at grid 2601. 
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Figure 13. Case 2, pitch acceleration time history at grid 2601. 
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Figure 14. Case 2, Z velocity time history at grid 2601. 
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Figure 15. Case 2, pitch rate time history at grid 2601. 
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Figure 16. Sensitivity study configurations for number of fibers. 
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Figure 17. Case 2, pitch rate time history at grid 2601. 


14 

American Institute of Aeronautics and Astronautics 



References 

Tgawa, H., Murayama, H., Nakamura, T., Yamaguchi, I., Kageyama, K. Uzawa, K., Wada, D., Ohsawa, I., Kanai, 
M., and Omichi, K., “Measurement of Distributed Strain and Load Identification Using 1500 mm Gauge Length FBG 
and Optical Frequency Domain Reflectometry,” Proceedings of 20 th International Conference on Optical Fibre 
Sensors , Vol. 7503, 2009. 

2 Cusano, A., Cutolo, A., Nasser, J., Giordano, M., and Calabro A., “Dynamic Strain Measurements by Fibre Bragg 
Grating Sensor,” Sensors and Actuators A: Physical , Vol. 1 10, No. 1-3, 1 February 2004, pp. 276-281. 

3 Ko, W. L., Richards, W. L., and Tran, V. T., “Displacement Theories for In-Flight Deformed Shape Predictions 
of Aerospace Structures,” NASA TP-2007-214612, 2007. 

4 Kang, L.-FF, Kim, D.-K, and Han, J.-J, “Estimation of Dynamic Structural Displacements Using Fiber Bragg 
Grating Strain Sensors,” Journal of Sound and Vibration , Vol. 305, No. 3, 2007, pp. 534-542. 

5 Moore, J. P., “Method and Apparatus for Shape and End Position Determination Using an Optical Fiber,” U.S. 
Patent No. 7,813,599, issued October 12, 2010. 

6 Shkarayev, S., Krashantisa, R., and Tessler, A., “An Inverse Interpolation Method Utilizing In-Flight Strain 
Measurements for Determining Loads and Structural Response of Aerospace Vehicles,” Proceedings of Third 
International Workshop on Structural Health Monitoring , 2001. 

7 Tessler, A., and Spangler, J., “A Variational Principle for Reconstruction of Elastic Deformations in Shear 
Deformable Plates and Shells,” NASA/TM-2003-2 12445, 2003. 

8 Tessler, A., Spangler, J., Mattone, M., Gherlone, M., and Di Sciuva, M., “Real-Time Characterization of 
Aerospace Structures Using Onboard Strain Measurement Technologies and Inverse Finite Element Method,” 
Proceedings of the 8th International Workshop on Structural Health Monitoring , Vol. 1, 201 1, pp. 981-989. 

9 Pak, C.-g., “Wing Shape Sensing from Measured Strain,” AIAA-2015-1427, 2015. 

10 O’ Callahan, J., Avitabile, P., and Riemer, R., “System Equivalent Reduction Expansion Process,” Proceedings 
of the 7 th International Modal Analysis Conference , Las Vegas, Nevada, 1989, pp. 29-37. 

n Pak, C.-g., and Friedmann, P. P., “New Time Domain Technique for Flutter Boundary Identification,” 
AIAA-92-2102, 1992. 

12 Doggett, R. V. Jr., Rainey, A. G., and Morgan, H. G., “An Experimental Investigation of Aerodynamic Effects 
of Airfoil Thickness on Transonic Flutter Characteristics,” NASA-TMX-79, 1959. 

U MSC/NASTRAN Quick Reference Guide Version 69, The MacNeal Schwendler Corporation, Newport Beach, 
California, 1996. 

14 Pak, C.-g., and Truong, S., “Creating a Test- Validated Finite-Element Model of the X-56A Aircraft Structure,” 
Journal of Aircraft, Vol. 52, No. 5, 2015, pp. 1644-1667. 

15 NASA Technical Standard, “Load Analysis of Spacecraft and Payloads,” NASA-STD-5002, June 21, 1996. 


15 

American Institute of Aeronautics and Astronautics 



